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ABSTRACT 

An analysis is made of current Navy procedures for 
assigning a probability density function to demand and 
the technique (s) used to forecast the parameter (s) of 
the particular density function chosen. The gamma den- 
sity function is investigated as a possible replacement 
for the three density functions currently being used. 

A comparison of the gamma and normal density functions 
is made with regards to inventory costs, observed protec- 
tion levels, and unit effectiveness. Additionally, a 
comparison is made, through a simple simulation model, 
of demand patterns generated by demand reporting as it 
is done today and consumption data reporting as it might 
be done in the future. The intent of the simulation is 
to provide some insight into the impact the demand re- 
porting method might have on the variance of demand. 
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INTRODUCTION 



i 

Today, the Naval Supply System faces the task of sup- 
porting very complex weapon systems under the burden of a 
funding climate that is already austere and will probably 
become more so. Accordingly, this situation makes effi- 
cient inventory management of utmost importance. This is 
not to say that inventory management has not been vital in 
the past, but rather the decision-making process with re- 
gards to inventory management becomes more critical. Ob- 
viously as the dollars the supply system receives are 
reduced, the application of the dollars must be judicious- 
ly applied to provide the best support possible to the 
fleet. By support it is meant the repair parts required 
to keep a weapon system or essential shipboard/aircraft 
operating. In view of the fact that an effective essen- 
tiality coding has not been established within the Navy, 
this means supporting almost everything with a Navy stock 
number. It should be noted that the general objective of 
the Naval Supply System Command for the past several years 
has been to reduce the time equipment is down or not 
operating due to the failure of a repair part. This of 
course entails more than having the part in inventory; 
however, this is the most critical aspect since most items 
required are unique to the Navy and are not readily avail- 
able from industry. 
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Essential to having the right part on hand at the right 
time is the ability to predict future demand for repair 
parts. This is by no means an easy task in view of the 
random nature of demand over a given time interval. The 
purpose of this thesis is to appraise what is being done 
within the Naval Supply System to predict future demand 
and to investigate possible alternative techniques. 

In predicting future demand there are two major factors 
that must be addressed. They are: 

1. Determine the probability distribution that repre- 
sents or at least approximates the actual distri- 
bution . 

2. Develop a technique for forecasting the parameters 
of the probability distribution determined to be 
representative of how future demand is expected to 
be generated. 

The current Navy methods of forecasting demand will be 
investigated. The selection of the probability distribu- 
tion will then be discussed with primary emphasis on com- 
paring the normal distributions presently in use with gamma 
distributions as to their effect on the establishment of 
inventory levels. A- brief discussion is made with regards 
to the use of the Poisson, negative binomial and the pseudo- 
normal distributions. The method of demand collection is 
then examined, and some possible problems associated with 
the current system within the Navy are discussed. Finally, 
the conclusions and recommendations generated by this study 
are presented. 
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II. DEMAND FORECASTING 



There are many ways to forecast demand; however, most 
of them can be grouped into two broad categories. The 
technical or subjective method and the scientific (some 
people would prefer mathematical) or the objective method. 
The technical method is the method by which the technical 
experts in a particular field give their best subjective 
estimate as to the demand for an item. This method appears 
appropriate only when there does not exist past demand 
data for an item. It could be useful possibly as a weight- 
ing factor in one of the scientific methods. Presently, 
the Navy supply system uses the technical method for com- 
puting Best Replacement Factors (BRF) for new items of 
supply or for items which have not realized any demand 
but are required to be stocked. This method will not be 
discussed any further, since the objective of this thesis 
is to deal with the items having recurring demand. 

The scientific method of forecasting demand breaks 
down to basically two types of technique. The first in- 
cludes the time series techniques where past demand pat- 
terns are used to predict future patterns. Some examples 
of this type are the straight average, moving average, 
single and double exponential smoothing, spectral analysis 
and Bayesian procedures. These techniques rely solely on 
looking at the pattern of past demand over time to predict 
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future demand for a specified time frame. The second type 
within the scientific method is the causal type. Although 
this type also relies on past demand patterns, it attempts 
to explain the causes for the particular patterns. There- 
fore, this technique requires the knowledge of certain 
causal variables which in turn can be related to the depen- 
dent variable for prediction purposes. 

This is the ideal technique; however, it is very time 
consuming, costly, and usually very difficult to determine 
the causal variables. This technique is usually related 
to econometric models. The Navy uses a regression model 
in determining the Mean Absolute Deviation (MAD) of demand, 
which is further used to determine the standard deviation 
of demand. In this model the explanatory variable or 
causal variable is the average demand D. This model will 
be discussed in more detail later in this section. 

The following time series techniques are defined: 

A. STRAIGHT AVERAGE — This technique does nothing 
more than divide the total number of demands by 
the number of observations as dictated by the 
desired time interval. For example, if average 
monthly demand was desired, the total number of 
demands would be divided by the total number of 
months for which demand data was available. 

B. MOVING AVERAGE — In laymen's terms this method 
drops the oldest observation of demand and adds 
the latest and then computes the straight average. 
The computational form selects a number of base 
observations, say four months, eight months, or 
four quarters. Once the base is selected, a 
forecast is generated by computing the average 

of the base observations. Then when a new 
observation is available, the new forecast is 
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computed by adding to the old forecast the new 
observation minus the oldest observation divided 
by the number of base observations. Mathemati- 
cally, 



NEW FORECAST 



OLD FORECAST + 



NEWEST OBS-OLDEST OBS 
Number of base obs . 



C. EXPONENTIAL SMOOTHING - This technique is used 

to assign different weights to different observa- 
tions of demand. In the straight and moving 
average methods equal weight was assigned to each 
and every observation. In exponential smoothing 
this is not the case. The more recent demands 
are given heavier weights than the older demands. 
Mathematically the forecast is calculated as: 

NEW FORECAST = OLD FORECAST + a (NEW OBS - OLD FORECAST) 
which can be further simplified to 

NEW FORECAST =(l-a)OLD FORECAST+ (a) NEW OBSERVATION 

where a = assigned weight to most recent observa- 
tion and 0 <_ a <_ 1. 

Now that the forecasting techniques have been defined 
it would be worthwhile to look at the advantages and dis- 
advantages of each technique. An advantage of the straight 
average is the fact that it provides the maximum likelihood 
estimate of the mean demand when a normal distribution de- 



scribes lead-time, demand. In other words, the sample mean 
is an unbiased efficient estimator of the actual mean of 
the distribution. The disadvantage of using the straight 
average is its slow reaction to large increases or decreases 
in demand. To use this latter point as an argument against 
the use of this technique seems to be contradictory to the 
assumption of lead-time demand being normal with a constant 
mean. Since this technique produces the best estimator of 
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the parameter, it would seem that the more observations 
available the closer this estimate would be to the actual 
parameter regardless of the immediate consequences of a 
significant change in a particular demand observation. 
Another disadvantage of the straight average lies in its 
inability to adjust to changes in demand patterns caused 
by a varying mean. 

The advantage of the moving average over the straight 
average is its ability to react more quickly to changes 
in demand patterns. The disadvantages associated with 
using this statistic are high computer costs in both space 
and time, which are due to the requirements to store all 
past observations over the base period. 

Exponential smoothing provides an estimate of the mean 
which is asymptotically unbiased and it reacts more quickly 
to changes in demand patterns. It also requires very 
little computer space and time. These latter attributes 
make it the most desirable technique as far as the Navy 
Supply System is concerned. 

Presently the Navy Uniform Inventory Control Proce- 
dures (UICP) uses exponential smoothing for forecasting 
demand; however, this forecast is refined by use of trend 
tests and a tracking signal filter. The trend test is 
used to verify if there is actually a change in demand 
pattern. The tracking signal is a statistic which is used 
to determine if the forecasting rule is performing well . 
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It is an algebraic sum of forecast errors. If the fore- 
casting rule is operating properly the sum of these errors 
should be close to zero; if it is large, then the system 
is considered out of control and recomputation is performed. 

So far, only techniques for forecasting the mean 
(average demand) have been discussed; however, the variance 
which is a measure of how demand varies about its average 
must also be forecasted. Estimating the variance accurate- 
ly is critical since it has a pronounced affect on the 
reorder point and safety levels which of course drive the 
stock levels maintained. The stock levels further determine 
the major portion of the inventory investment cost. Pres- 
ently the UICP model of the Navy estimates the variance by 
using the mean absolute deviation (MAD) . The MAD is defined 
as the expected value of | X— n | and the most commonly used 
estimate of MAD is the sample average of the absolute differ- 
ences between the observed demands and their sample mean, 

N |x.-x| 

Z 



i=l 



N 



The UICP model does not use this formulation but rather 

-'N _g 

uses the formula, MAD = aD p , where D represents expected 
lead-time demand. The formula is based on a study done by 
Fleet Material Support Office (FMSO) where regression analy- 
sis was used to determine the coefficients a and 3* This 
study [1] was done in 1962, and the values of a and 3 were 
determined to be 1.37 and .717, respectively. 
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The UICF model takes the MAD computed as above and 
then uses as the estimate of the standard deviation 

A 

S = 1.25 MAD. Once a MAD is determined for an item 
successive MADs are exponentially smoothed the same as 
the estimate of mean demand. The use of this technique 
for estimating the standard deviation is based on work 
done by R. G. Brown [2] . This technique has come under 
criticism by several investigators including particularly 
Dr. P. W. Zehna [3] when used with a probability density 
function other than the normal. 

Using a 1000 item sample of demand data provided by 

the Air Force, which has been used in studies done by all 

services in the past several years, a 50 item subsample 

was selected randomly and estimates of the standard devia— 

E I x— x I 

tion and MAD (MAD^ = — 1 — ^ — L form) were computed. The 
sample standard deviation was then compared with the value 
1.25 (MAD^) and a sign test performed to test the hypo- 
thesis that there was no difference between the two esti- 
mates of standard deviation. Since the difference between 
the two standard deviations was always positive, the sign 
test strongly rejected the hypothesis. 

In view of the fact that the translation between the 
theoretical MAD and the standard deviation appears to be 
questionable it seems prudent to also investigate the use 
of the linear regression technique of estimating the MAD. 
Problems encountered with this technique would certainly 
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compound any problem associated with utilizing the MAD to 
estimate the standard deviation as is presently being done. 

Great care must be taken when using a linear regression 
prediction scheme where the values of the explanatory 
variable (the predictions of mean monthly demand) used in 
estimating the least squares coefficient are not known. 

Using a regression model to predict beyond the range of 
values of the explanatory variable can lead to serious dif- 
ficulties. From the random sample of 50 items, the regres— 

sion estimate of the mean absolute deviation, MAD = 1.37 
— 717 

x' , was calculated for each item. These values were 
compared with the other estimate of the mean absolute devia— 

tion, MAD^ . It was noticed that MAD was always greater than 

/\ 

MAD^ when MAD^ was less than one and MAD tended to be less 
than MAD^ when the latter value exceeded one. The hypothesis 
that the two estimates of mean absolute deviation are the 
same was tested using the Wilcoxon Signed Rank Test, and 
it was strongly rejected at the 0.01 level of significance. 
Results of that test are found in Appendix A. 

The above results combined with the earlier comparison 
of the sample standard deviation and 1.25 MAD-^ point out 
the possible source of problems experienced by the Navy in 
estimating variance. The wide disparity between the differ- 
ent estimates of the standard deviation can best be illus- 
trated by presenting the three estimates for each of the 
50 items. This is done in Table I. 
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TABLE I 



Comparison of three methods for estimating 
standard deviation of demand. 
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1.25 


0.26 
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15.48 


12.25 


C . 72 


0.54 


544.56 


247.55 


C . 12 


0.04 


43.53 


29.95 


16.20 


12.08 


2.58 


2.43 


2.55 


1.52 


6.48 


4.54 


C . 54 
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26 C .23 


163.49 


3.46 


1. 77 


57. 17 


35.46 


11.52 


5.54 


0.82 


0.37 


1 . 16 


C .65 


2.23 


1.70 


1 . 16 


1.02 


3.50 


2.31 


1.14 


C . 77 


1 . 17 


C . 73 


0.81 


0.45 


107 . 11 


. 50.10 


6.48 


3.70 


16.40 


7.53 


4.27 


2.16 


C . 5 1 


C . 66 


C . 52 


0.32 


29.50 


21.41 


3.79 


3.57 


4. 05 


3 .35 


0.26 


0.16 


4.65 


3.60 


3.83 


2.72 


27.55 


17.48 


4.71 


3.93 


5.81 


5.54 


26. 12 


17.08 


4 .47 


3 .67 


13.22 


6.22 


18.88 


7.03 


32.72 


22.08 


0.85 


0.42 


1.21 


1.06 


166.55 


90.25 



MAD 1 1.25 MAD 

0.26 
3 .58 
0.26 
0.58 
0 . 16 
11.69 
0.66 
57 .23 
0. 1C 
13.89 
6.28 
2 . 14 
1.35 
3 .15 
0.65 
53.61 

1 .50 
25 . C 5 

3 .60 
0.46 

0 . 65 
1.53 
1.24 

1 . 75 
0.81 
0 .75 
0 .52 

47 . 12 
2.62 

4.71 

2 .32 
0.79 
0.42 
9.56 
3 .00 
2 ,< j 4 
0.26 

2.72 
2.43 
7.55 
3 . 5 C 
5 .50 
7.87 
2 .96 
3 . 63 
3.5 1 

10.05 
0 . 5 C 
1 .14 
23.77 
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The method of forecasting the two critical parameters, 
the mean lead-time demand and the standard deviation of 
lead-time demand, discussed in this section will also be 
compared in Section III in order to appraise the effects 
on inventory levels caused by their use. 
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III. SELECTION OF PROBABILITY DENSITY 
FUNCTIONS TO REPRESENT LEAD-TIME 



DEMAND 



Whereas there has been a great deal of concern in the 
past several years regarding the ability to forecast lead- 
time demand and variance, little or no concern has been 
evidenced regarding the demand distribution itself. Ob- 
viously no one probability distribution function is going 
to fit the demand for all military items; however, one is 
interested in determining if there is any one distribution 
that can best represent demand as the real world sees it. 

An affirmative answer to this question would require a ver- 
satile distribution function since it is readily recognized 
that demand for military items covers the range from very 
low to very high, and frequently is erratic in nature. The 
second question to ask, given the first is answered affir- 
matively, would be the feasibility of using such a distri- 
bution. Three probability distributions are used in UICP 
in varying degrees to describe demand. These distributions 
are the Poisson for very low demand items (i.e. average 
annual demand £ 1) , the negative binomial for items with 
annual demand in the range of two to ten, and normal for 
items with average annual demand greater than ten. 

Prichard and Eagle [4] state that empirical demand dis- 
tributions for items with low means are usually skewed to 
the right thereby making is desirable to use the Poisson 
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and negative binomial distributions to describe this type 
of item. They further state that the Poisson, negative 
binomial, and normal distributions seem to satisfactorily 
approximate demand for a majority of items and also indi- 
cate that price and consumer behavior usually are the 
causes for differentiating between the use of the Poisson 
and the negative binomial. That is, a cheap low demand 
item will be ordered in batches thereby causing a high 
variance and suggesting the use of the negative binomial. 

On the other hand, an expensive low demand item will be 
ordered one at a time thereby being better represented by 
the Poisson. Although not specifically stated, the reader 
is led to the assumption that the high demand items can be 
represented by the normal regardless of the price. It is 
by this reasoning that they contend a single mathematical 
distribution cannot represent all demand. They do not 
address, however, the possibility of a single distribution 
satisfactorily approximating these three distributions. 

As mentioned previously, the Navy Supply System is using 
these distributions in varying degrees. A recent study 
completed by FMSO indicates that one Navy Inventory Control 
Point (ICP) is using the normal to represent all of its 
demand [5] . Additionally, the use of the Poisson has de- 
creased in the past several years, as tables have been 
generated to use in its place. Basically when discussing 
demand distributions within the Navy it is the normal 
distribution that is being referred to. 
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Recently there has been a move underway to substitute 
the pseudo— normal distribution for the normal and the 
other distributions. The pseudo— normal has been recom- 
mended because of the computational advantages that it 
allows in the calculation of inventory reorder levels and 
reorder points. Because it performs about the same as 
the normal with regard to describing demand patterns, we 
will concentrate our attention on the normal distribution 
in this study. 

As stated previously, representing military demand by 
a single probability distribution would require a very 
versatile distribution. It is felt that the gamma family 
of distributions best meets this qualification. Dr. 

Peter W. Zehna [6] states: "The family (gamma) is so ex- 
tensive in shapes of densities available that it is a 
fairly safe assumption to make as a model for an experi- 
ment described by almost any non— negative random variable." 
Additionally Hoel, Port, and Stone [7] point out that in 
most cases involving a random variable X which is known 
to be positive the assumption that X has a gamma density 
will provide an approximation or at least an insight into 
the true but unknown situation. Finally, the versatility 
of the gamma is shown graphically in Figure 1 where curves 
associated with the Poisson, negative binomial, and normal 
distributions are displayed in 1A and the curves in IB 
represent a particular set of gamma distributions obtained 
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by varying the parameters A and a. These parameters re- 
ferred to are called the shape (A) and scale (a) para- 
meters of the gamma function and will be discussed in 
more detail later. 




Figure 1 

A true but unknown situation seems to describe more 
than adequately the situation inventory managers find 
themselves in when trying to predict how demand is going 
to occur. Therefore it seems appropriate that the gamma 
should be investigated as to its ability to represent 
demand of military items and to be compared with the normal 
as to which better describes the situation. It is pre- 
cisely this objective that the remainder of this section 
is devoted to. 

Using the sample of 50 items chosen from the 1000 item 
sample of Air Force data, the monthly mean (X) , standard 

A 

deviation (S) and mean absolute deviation (MAD) were 
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computed based on 57 months of past demand. Additionally, 
the gamma parameters A and a were computed using the 
method of moments in the following manner: 

FOR T(x,A,a)= (ax) A ^e aX ; E(x)= A/a, V(x)= A/a 2 

THEN solving for A and a 

x/S 2 = a and x 2 /S 2 = A 

These parameters were computed mechanically by the computer 
which at the same time was performing a goodness of fit 
test for two different cases. In the first case the normal 
distribution was fit to the data using the sample mean and 
the sample standard deviation which are the maximum likeli- 
hood estimators for the parameters of the normal distribu- 
tion. In the second case the gamma distribution was fit 
to the data using the sample mean and variance to compute 
the estimates of the gamma parameters A and a. A Kolmogorov— 
Smirnov (K— S) test was performed in each case resulting in 
a strong rejection of the normal hypothesis and an accept- 
ance. of the majority of items under the gamma hypothesis. 
Table II provides a summary of the results of the K— S test. 

TABLE II 

Results of Kolmogorov— Smirnov test on 
50 randomly selected items. 



Distribution 


No. Items 


No. Items 


Percent of Items 


Fitted 


Accepted 


Rejected 


Accepted 


— ? 

Normal (x,S ) 


0 


50 


0% 


Gamma (A, a) 


37 


13 


74% 
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A listing of results obtained in both cases is provided 
in Appendix B. 

The necessity of fitting a probability distribution 
to demand data comes from the requirement to protect 
against running short of stock. In inventory models pro- 
tecting against shortages is achieved by ordering stock 
before the amount on hand falls too low. The procurement 
policies specify that orders will be placed as soon as 
the stock assets reach or fall below a given level called 
the reorder level. The higher the reorder level, r, the 
lower will be the probability that the total demand dur- 
ing the lead-time exceeds r. Define the risk to be this 
probability and let its complement be called the protection 
level. 

RISK = Pr [Total lead-time demand > r] 

The reorder level is adjusted up or down with the degree 
of adjustment depending on the desired protection level 
and the probability distribution of lead-time demand. The 
amount of material ordered, the reorder quantity, is usually 
constrained by the administrative cost to place an order, 
the cost of holding that amount of material in inventory 
and the cost to the system of being out of stock when a 
demand occurs. These costs are normally referred to as 
ordering, holding, and shortage costs and are not easily 
obtainable. Consequently they are usually a product of 
the particular inventory system. 
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The Navy uses a modification of Hadley and Whitin' s lot 
size— reorder point model (<Q,r> model) defined as follows: 



1. Let I be the annual holding cost rate per unit of 
stock, Ax the fixed cost of placing an order, C 
the cost per unit for item i, Tj_ mean annual demand 
for item i, and it the shortage cost for each unit 
short. (A fixed shortage cost is not assigned for 
each item; rather, the penalty ir is simply mani- 
pulated to meet budget and risk constraints.) If 
H(r) is the probability that the total lead-time 
demand for item i will exceed r, the reorder level 
is determined from the equation 



H(r i ) 



IC.Q 

I 

IC . Q . + ITT 

II 



2. The reorder quantity, Qi , for item i is taken to 
be the maximum value among (1, Qi , x^/4). Where 




2x^Ax 

IC. 

l 



is the well known Wilson lot size. 

It is important to note that under the Navy UICP model the 
computation of is completely independent of the reorder 
point and the risk, thereby making independent of the 
demand probability distribution. 

Prior to comparing the UICP model under the assump- 
tions of normal versus gamma distribution of lead-time 
demand it is necessary to obtain a feeling for the rela- 
tive difference between the theoretical reorder points 
assigned by the gamma versus those assigned by the normal. 
For an illustrative purpose three items were chosen from 
the 50 item sample. Table III shows the results of 
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comparing the empirical, gamma, and normal distributions. 
Appendix D gives the graphic representations of the 
empirical cumulative distributions from which the actual 
protections were determined. The two normal distribu- 
tions differ only in the choice of standard deviations. 

In one case the sample standard deviation is used and in 
the other case SD is computed using the regression for- 
mula for the MAD as discussed in Section II. The two 
sets of figures for each distribution represent 80 and 50 
percent protection levels respectively. Table III indi- 
cates the gamma does a much better job at approximating 
the empirical distribution for the first two items than 
does either of the normal distributions. For the third 
item the gamma, and the normal using SD, both do a good 
job of approximating the empirical distribution at the 80 
percent protection level. At the 50 percent protection 
level only the normal with SD does well. At this level 
the gamma underprotects as much as the normal using S over- 
protects. Looking at a histogram of this item (Figure 2) 
it can be seen that the gamma is not sensitive enough to 
the three large monthly demands (662, 757, 1733) and the 
normal is too sensitive. This situation seems to be the 
crux of the inventory problem. Should the system over- 
react to peaks in demand or should it ignore them? Ob- 
viously it costs inventory dollars to over-react such as 
in this case of stocking 71 more units than required; 
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S = Sample standard deviation 

SD = Computed standard deviation (SD = 1.25*MAD) 
— 717 

MAD = 1.37 x* /X/ 



Histogram of Item #16 (AF #651) 
x = 120.5 S = 258.0 
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Median Gamma Median Empirical „ Mean > Median, Normal 

r Figure 2 



however, how does one measure the cost of the shortages 
one can anticipate by only stocking 16 units. It is in- 
teresting to note the effect the computed standard 
deviation, SD, has on this item. The value for SD is 
approximately 20 percent of the value S. This is not 
difficult to understand when one considers the fact that 
in the computation of the MAD the mean is raised to the 
.717 power and then multiplied by 1.37. Except for mean 
monthly demands less than 3.04 this formulation will al- 
ways give a standard deviation less than the mean. 

The final comparison between the gamma and the normal 
is made with regards to the inventory levels and the 
associated costs that are generated under each distribu- 
tion. A FORTRAN program was written to compute reorder 
levels and reorder quantities for each of the 50 items 
chosen from the demand data base. These values were then 
used to generate holding costs. In addition, shortages 
were accumulated and an overall unit effectiveness was 
computed. Overall unit effectiveness is defined as the 
total number of units satisfied divided by the total num- 
ber of units demanded. Appendix C contains three programs 
to accommodate the computation of the gamma cdf, and the 
use of the exponentially smoothed average and the predicted 
standard deviation for the normal. Two different runs 
were made to determine the impact that different parameters 
would have on the levels computation and associated costs. 



27 



The procedures of the program for the first run are as 
follows : 

1) Read the following values: 

Holding Rate = XI = . 01/month 
Ordering Cost= AX = $25. 

2) Generate unit price by randomly selecting values 
between one cent and 50 dollars. The same seed 
for the random number generator was used to 
obtain the same sequence of unit costs each time 
the program is run. 

3) Compute the means, variances, and standard devia- 
tions using the first 21 months of demand history. 

4) Compute the initial reorder levels and quantities 
and randomly set the initial assets on hand between 
r+1 and r+Q. It is in this step where the com- 
putation of r is affected by both the probability 
distribution and the parameters chosen and the Q 

by the particular forecasting technique chosen 
(i.e. sample mean, exponential smoothing, etc.). 

5) Run the program for the remaining 36 months compar- 
ing each month's actual demand with the available 
stock on hand and computing units short, if any, 
total holding cost and finally overall unit effec- 
tiveness. 

6) The final step is to compute the observed protec- 
tion. This is done by dividing the total number 
of lead-time periods in which shorts do not occur 
by the total number of lead-time periods. 

The first case was run three times for protection levels of 
30, 50, 90 percent for each of the three distributions; 
the normal using the sample mean and sample standard devia- 
tion; the normal using the exponentially smoothed average 
and the exponentially smoothed estimate of the standard 
deviation; and finally the gamma using the parameters esti- 
mated by the method of moments. The results are provided 
in Table IV. 
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TABLE IV 



OBSERVED PROTECTION, TOTAL HOLDING COST, TOTAL UNITS SHORT, 
AND OVERALL UNIT EFFECTIVENESS FOR VARIABLE r and Q - MAX 
(1.0, Q , t/4 ) . 



Desired Protection Level = 30% 





Normal 

(x,S) 


Gamma 
(A, a) 


Normal 
(SMO , SD) 


Observed 

Protection 


25.81% 


19.67% 


39.29% 


Total Hold- 
ing Cost 


42,379 


42,217 


55,734 


Total Units 
Short 


1,661 


1,732 


1,311 


Overall Unit 
Effectiveness 


85. 73% 


85.12% 


88.74% 



Desired Protection Level = 50% 



Observed 

Protection 


43.75% 


31.15% 


56.90% 


Total Hold- 
ing Cost 


57,075 


42,897 


61,236 


Total Units 
Short 


763 


1,474 


845 


Overall Unit 
Effectiveness 


93.44% 


87.33% 


92.74% 



Desired Protection Level = 90% 



Observed 

Protection 


59.38% 


58.46% 


70.00% 


Total Hold- 
ing Cost 


107,195 


67,442 


74,709 


Total Units 
Short 


454 


498 


409 


Overall Unit 
Effectiveness 


96.10% 


95.72% 


96.49% 
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Significant among the results depicted in Table IV are 
the high unit effectiveness corresponding to the low ob- 
served protection levels. In fact the observed protections 
fell short of the desired levels in all but two cases where 
the normal distribution was used with the parameter esti- 
mates SMO and SD at the 30 and 50 percent desired protec- 
tion levels. These results point out that the system can 
satisfy most of the demands (high unit effectiveness) while 
achieving very low protection levels. Both the short lead- 
time (taken to be one month in this evaluation) and the 
large values of the reorder quantities Q relative to the 
reorder levels r account for this phenomenon. For in these 
circumstances the system rarely relies on the reorder level 
to protect against stockouts. In most cases the available 
stock is greater than the reorder level. 

Both the protection level and the unit effectiveness 
figures obtained by the gamma distribution fall below the 
figures obtained by the normal distributions. At first 
this seems an indictment against use of gamma distribution. 
However, one must consider the extra cost of achieving this 
higher protection. That is, one must ask if the trade off 
is cost effective. 

To obtain higher protection levels the normal distri- 
butions must be choosing higher reorder levels on the 
average. These higher protection levels result in higher 
investment and holding costs. Thus a comparison of the 



holding costs gives an indication of the cost of achieving 
the higher protection. Although no definite conclusions 
can be drawn from the values in Table IV there are indica- 
tions that the gamma performs better than the normal with 
the same estimates of the mean and standard deviation at 
the higher protection levels when measured in terms of unit 
effectiveness per dollar of holding cost. Also, when com- 
paring the two normal cases, the one with the exponentially 
smoothed estimate appears to perform better than the other 
using unit effectiveness per dollar of holding cost. 

In order to better assess the effects of the type of 
probability distribution on the measures of effectiveness 
and the costs, the reorder quantities were all forced to 
be Q=l. This allows us to focus more clearly on the impact 
of the reorder levels determined by the two types of dis- 
tribution. Also only the two cases where the sample mean 
and the sample standard deviation were used to estimate 
the parameters were considered so that any differences 
could be contributed to the fits of the probability distri- 
butions. Table V shows the results obtained with the re- 
order quantities set at one. Only the observed protection 
levels and the holding costs are presented in Table V 
because when Q=l, the protection level should be a better 
measure of performance than unit effectiveness. 

The dominance of the gamma over the normal is better 
illustrated by Figure 3 . There it is clearly shown that 
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TABLE V 



OBSERVED PROTECTION AND HOLDING COSTS 
AT VARIOUS LEVELS OF PROTECTION. 



Desired Protection Level 


30% 


Demand Distribution 


Normal (x,S) 


Gamma (A, a) 


Observed Protection (OP) 


25.6% 


23.5% 


Total Holding Cost (THC) 


1,541 


1,399 




C 


>0% 


OP 


62 . 5% 


33.6% 


THC 


11,658 


2,265 




, 60% 


OP 


68.7% 


42.9% 


THC 


19,868 


3 ,352 




70% 


OP 


75.5% 


56.2% 


THC 


29,018 


5,876 




80% 


OP 


79.2% 


69.3% 


THC 


40,017 


12,377 




90% 


OP 


83.2% 


81.0% 


THC 


55,480 


27,154 
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Figure 3. Holding Cost (,000) 



the gamma protection is greater than the normal protection 
for any value of holding cost. These results coupled with 
the comparisons involving unit effectiveness versus dollars 
of holding cost support the hypothesis that the gamma is 
more cost effective than is the normal. 

The values obtained for actual protection in Tables IV 
and V appear to be inconsistent since one would expect 
greater protection with higher inventory levels. The prob- 
lem appears to be in the definition of protection and 
relates to the frequency at which the system reaches its 
reorder point. The effect that Q can have on the protec- 
tion level is best shown by an example. A year's demand 
for one of the items of the 50 item sample is recorded in 
Table VI. For ease of computation the variable Q was com- 
puted as the total yearly demand divided by 4(Q=x/4) and 
the reorder level was calculated from a normal distribution 
to give a 50 percent protection level. The initial on hand 
level of stock was assumed to be r+1. 



TABLE VI 

PERFORMANCE COMPARISON OF Q=1 vs. Q=x/4 . 



Demands 


4 8 3 0 2 14 117 43 7500 


x = 203 


Q=x/4=51 




x= 17 


OH 


14* 57 54 54 52 38 -79* -20* 24 19 19 19 
Protection = 1 — 2/3 = 33% 

U.E. = 1 - 99/203 = 55% 


r= 17 
OH=18 


Q = 1.0 






OH 


14* 10*15*18 16* 4*— 94* -25* 11*13*18 18 
Protection = 1 — 2/9 = 77% 

U.E. = 1 - 124/203 = 39% 





where * implies an order was placed; U.E.=unit effectiveness 
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Finally, a discussion of the overall unit effective- 
ness figure used in Table IV is warranted. An important 
question being asked today is whether an inventory system 
should be geared to protect against total units short or 
the total number of requisitions short. Obviously, if 
only one unit were demanded per requisition (which is an 
assumption of the Q, r model) these two measures would be 
the same. However, when examining demand patterns for 
the items used in this study it was noticed that there 
would be many months of demand within a particular range 
and then two or three months of very high demand. For 
example, demand in the range of 50 to 100 would be received 
for several months and then a demand for 3300 would be re- 
ceived. It would be difficult to accept the assumption 
that the 3300 units demanded represented 3300 requisitions. 
The importance of this point is that if the 3300 units 
demanded were not on separate requisitions the value of 
unit effectiveness becomes distorted. 
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IV. METHOD OF RECORDING DEMAND 



Up to this point the problems associated with fore- 
casting demand have been defined as the ability to describe 
the probability density function of demand and the selec- 
tion of the proper forecasting scheme to predict the para- 
meters of that density function. One of the major factors 
contributing to the complexity of these problems is the 
nature of demand patterns. An examination of the 1000 item 
sample used in this study indicated that in the majority of 
cases the demand patterns were very erratic. The primary 
effect of erratic demand is to cause large variances thereby 
rendering the task of fitting a probability distribution to 
the demand data, particularly difficult. More important to 
the inventory manager, however, is the effect these large 
variances have on his inventory levels. 

One of the major factors causing demand patterns as the 
inventory manager- sees them is the demand reporting system 
that., is utilized. If a ship requires a repair part, then 
the demand is created at the shipboard level at a particular 
point in time. Under current procedures there are four 
echelons of supply from which the ship can obtain the repair 
part required. Therefore the initial point at which the 
demand is recorded depends on where the requirement is 
satisfied. The four echelons referred to are the ship 
itself, a mobile logistics support ship, a stock point and 
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finally, an inventory control point. These echelons are 
critical in determining when a demand gets recorded at 
the ICP, which is the inventory manager and has the res- 
ponsibility of setting system inventory levels. For 
example, if a shipboard demand occurs on the first of the 
month and the ship can satisfy the demand from its own 
stocks then the demand is recorded at the shipboard level 
only at that point in time. 

If the demand did not reduce the onboard stocks to the 
ship's reorder point this particular demand would not be 
recorded at the stock point until sufficient demands 
occurred to reduce the stock to the reorder point, at which 
time the ship would place a replenishment demand on the 
stock point. Another factor affecting the time at which 
the demand is reorded at the stock point is whether the 
ship is at sea or in port since the ship has no means of 
communicating (except for high priority requirements) 
directly with the stock point if it is at sea. The point 
here is that the demand that a stock point observes could 
represent demand for that day or an accumulation of demand 
over a period of time. Since stock points report daily 
to the ICP's this is also the situation that the ICP in- 
ventory manager faces. There is such a large population 
of ships that enter and leave port at different times it 
would seem that over a period of time, the system demand 
would smooth itself. However, a review of actual data 
patterns does not- substantiate this. Thus, an investigation 
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into the method of collecting data seems warranted. Since 
time would not permit a thorough investigation it was de- 
cided to look at the effects of this situation through a 
simple simulation model in which demands for a single item 
were generated from seven different ships. It was assumed 
that the time between demands is exponentially distributed. 
The time at sea was assumed to be uniformly distributed 
between 0 and 45 days. Based on these assumptions the days 
on which demands occurred were randomly generated for each 
ship. Also the dates of arrival and departure from port 
were randomly generated, and initially all ships were con- 
sidered at sea. It was also assumed that all demands in- 
curred while the ship was at sea are reported to the stock 
point upon arrival in port. Additionally, any demands 
occurring while the ship was in port were reported to the 
stock point on the day they occurred. 

The alternative method of recording demand investigated 
in this simulation had the ship report the demand directly 
to the ICP on the' day it occurred. Inherent in this alter- 
native is the assumption that the ship has the capability 
of reporting demand directly to the ICP which is recognized 
as unrealistic at this time. However, since the ships 
currently keep records of daily demand through the 3M 
(Material, Management and Maintenance) system, it was felt 
that the above alternative was possible. Thus, we are com- 
paring consumption data vice demand data as currently 
reported in the Navy. 
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The simulation was run for a two year period with the 
demand recorded by month at both the ICP and stock point. 
Additionally, the mean and variance of demand over this 
two year period were computed for each reporting system. 
The results of this simulation are presented in Table VII. 

TABLE VII 

SIMULATED MONTHLY DEMAND 
AS RECORDED AT ICP AND NSC 



DATA 


YEAR 








MONTHLY 


DEMAND 










Consumption 

Data 


1 


9 


18 


18 


15 


17 


17 


18 


19 


15 


15 


23 


22 




2 


19 


11 


16 


12 


17 


18 


4 


12 


20 


12 


15 


8 


Demand 

Data 


1 


4 


7 


30 


4 


11 


16 


23 


23 


20 


17 


20 


24 




2 


7 


17 


9 


13 


8 


21 


2 


4 


40 


13 


22 


8 



Mean of Consumption Data = 15 Mean of Demand Data = 15 

Variance of Consumption Variance of Demand 

Data = 32 Data = 84 



As expected, the alternative of reporting demand direct- 
ly to the ICP (consumption data) shows a smoother pattern 
than that reported to the stock point (demand data) . The 
critical result however is that the variance of demand at 
the stock point was approximately times greater than the 
variance of demand recorded at the ICP. With smoother or 
more regular demand patterns and smaller variances the task 
of forecasting future demands should be easier. It would 
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be presumptuous to draw any conclusions based on such 
tentative results. However, further investigation into 
the use of consumption data does seem appropriate. 



40 



V. CONCLUSIONS 



A. DEMAND FORECASTING TECHNIQUES 

The following recomendations are made with regard to 
the forecasting schemes to be utilized to predict the 

mean and variance of the lead-time demand. 

2 

1) That the sample variance (S ) be used for predic- 
ting variance of lead-time demand regardless of 
the probability distribution chosen. There 
appears no statistical justification for using 
the Mean Absolute Deviation regardless of the 
particular estimate of MAD used. Since there 
was an apparent strong relationship between the 
sample MAD and the mean quarterly demand when 
the study was conducted at FMSO in 1962 it might 
be worthwhile to perform such a study with cur- 
rent demand. The use of the coefficients of 1.37 
and .717 computed in 1962 certainly does not 
appear warranted. 

2) The results of this study show no distinct pre- 
ference for the use of a sample average or the 
exponentially smoothed estimate of the mean. 



B. PROBABILITY DENSITY FUNCTION 

It is recommended that the current procedures of using 
three probability density functions to represent lead-time 
demand be replaced by the procedure of using the single 
distribution of gamma. This recommendation is based on the 
following conclusions: 

1) The gamma exhibits the versatility to represent 
low demand items as well as high demand items. 
Additionally the gamma is a non— negative probability 
density function which ensures that all of the 
probability mass will be represented by actual 
demand (i.e. no such thing as a negative demand 
such as is the case with the normal) . 
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2) The gamma appears to be less sensitive to peaks in 
demand and therefore consistently assigns lower 
reorder levels creating much lower inventory hold- 
ing costs. The argument here is that unless the 
demand pattern is smooth, stocking material under 
the normal assumption tends to put higher reorder 
levels than required to meet most demands. These 
higher reorder levels naturally result in higher 
protection levels but at the same time they require 
the system to carry a lot of stock and therefore 
incur high holding and investment costs to protect 
against infrequent demands. 

3) Finally, the conversion from the present system to 
the gamma would require little effort since the 
incomplete gamma function which gives the comple- 
mentary cumulative distribution is already pro- 
grammed and tabulated. Additionally using the 
method of moments to compute the gamma parameters 
requires minimum machine effort. 
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APPENDIX A 



VALUES OF S, SA, MAD^ AND MAD 
USED IN SIGN, AND SIGNED RANK TEST 



s 


SA 


C • 2 6 


0. 16 


7.56 


5.71 


C . r 7 


0. 17 


0.80 


C. 78 


C . 19 


o. ce 


15.48 


13.25 


C . { 2 


0 . 54 


544 .56 


247.55 


C .12 


0.C4 


4 2.52 


29.55 


1 6 . j C 


12. C8 


2.58 


2.43 


2.55 


1 . 52 


6.^8 


4.54 


C .54 


0 . c C 


26C .23 


163 .49 


2 . 46 


1 . 77 


57.]7 


35.46 


11.52 


5.54 


C * 


0.27 


1.16 


0.65 


2 . c 3 


1.70 


1 . 16 


1 . C2 


2.50 


2.31 


1. 14 


0.77 


1 . 17 


0.73 


C • £ 1 


0.45 


107.11 


90. 1C 


6.48 


3.70 


16. 4C 


7.53 


4.27 


3. 16 


C . 5 1 


0 . 66 


C . 5 2 


0.22 


29.5 0 


21.41 


2.79 


3.57 


4 . C5 


3.25 


C . 26 


0. 16 


4 . ( 5 


3. 60 


3.82 


2 . 72 


27.15 


17.48 


4.71 


3.93 


5.81 


5 . 54 


26. 12 


17 . ce 


4.47 


3.67 


1 3 • ^ k 


6 . 22 


16.88 


7. 02 


3 2.72 


22. CC 


C .65 


0.42 


I • 1 


1 . C6 


166.55 


90.25 



MAD X 


MAD 


0 . 13 


C . 2 1 


4.57 


3.19 


0 . 14 


0.21 


0.62 


0.79 


0.07 


0.13 


10.60 


9.35 


0.43 


0.53 


198. C4 


45.86 


0 . 03 


0.08 


23.96 


11.11 


9.67 


5.02 


1.95 


1.7 1 


1.22 


1.11 


3.63 


2.52 


0.48 


0.56 


130.79 


42.89 


1.42 


1.20 


28.37 


20.07 


4.43 


2.88 


0.29 


0.37 


0.52 


0.56 


1.36 


1.22 


0.82 


0.99 


1.85 


1.43 


0.62 


0.65 


0.58 


0.63 


0.36 


C . 42 


72.08 


37.70 


2.56 


2.10 


6.02 


2.77 


2.52 


1.86 


C .53 


0.63 


0.26 


0.34 


17 . 12 


7.97 


2.85 


2.40 


2.68 


2.11 


0 . 13 


0.21 


2 .88 


2.18 


2.18 


1.95 


13.99 


6.39 


3 . 15 


2.80 


4.43 


4.40 


13.66 


6.30 


2.54 


2.37 


4.58 


2.06 


5.63 


3.13 


17.66 


8.04 


0.33 


C .40 


0.85 


C .9 1 


72 .20 


19.02 



Results: Sign Test -Ho: S = SA. Since (S-SA) > 0 in all 

cases strongly reject hypothesis. 

Wilcoxon Signed Rank Test - Ho: MAD, = MAD 
Test statistic is T = —171.5 for negative sign at a = .01 
T c = 612.5 from Ref. [8] . Therefore strongly reject 
hypothesis . 
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APPENDIX B 



KOLMOGOROV-SMIRNOV TEST 



NORMAL 




GAMMA 


PROB = 0% 


PROB =0% 




PROB = 99.6% 


PROB = 100% 


0 


0 




100 


100 


0 


3.2 




100 


100 


0 


0 




6.6* 


62.6 


0 


0 




100 


100 


5.8 


0 




44.9* 


1.9' 


0 


0 




65.2 


99 


0 


0 




0* 


45.8' 


0 


0 




100 


100 


0 


0 




96.2 


99.7 


0 


0 




66.6 


88 


0 


0 




99.5 


99.9 


0 


0 




100 


99.6 


0 


0 




100 


55* 


0 


0 




98 


69.8 


0 


0 




1.8* 


92.6 


0 


. 8 




100 


72.7 


0 


0 




.5* 


42.2 


0 


0 




18.1* 


100 


0 


0 




100 


16.7' 


0 


0 




100 


100 


0 


0 




100 


95.5 


0 


0 




122* 


100 


0 


0 




100 


72.9 


0 


0 




99.7 


17.3' 



PROB = Probability of being incorrect if hypothesis is 
rejected 

* indicates gamma hypothesis is rejected 
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APPENDIX C 



FORTRAN PROGRAMS ASSOCIATED WITH COMPUTATION OF INVENTORY 
COSTS AND UNIT EFFECTIVENESS 

FORTRAN programs used in calculating values for Tables 
IV and V. 

The three programs in this Appendix represent computa- 
tions using a variable Q. For the values on Table V the 
programs were modified to set Q = 1.0. Additionally these 
basic programs were run with a counter to compute the 
number of periods when shortages occur and the total number 
of periods when orders were placed. This later modifica- 
tion was necessary to calculate the observed protection 
levels. 

Important labels used in the programs and not previously 
defined are as follows: 

BAR Mean monthly demand 

VAR Variance of monthly demand 

PT Desired protection level 

TSHT Total units short 

THCT Total holding cost 

OUE Overall unit effectiveness 



45 



FORTRAN PROGRAM FOR GAMMA (A, a) 
R E AL*8 A, XX 

DIN-ENSIGN IDEM (57 ) , DEM (57) , Z(30C ) 

A X = 2 5 . 

XSUM=0 . 

X SET = 0 . 

T HCT = 0 . 

X 1= .0 1 

20 R EAD( 5 * 10 , END=50 ) ( I DEM (I) , I = 1 , 57 ) 

IC FORMAT ( 121 5) 

DC 500 1=1,57 
5 C C CEN(I) = IDEM( I) 

P T= . 5 
Z ( 1 ) = 0. 

N = 2 1 
M =N 
K = 21 
H C T=0 » 

T SPT = 0 . 

T SUM=0 . 

DEV=0, 

I X= 1237 56 1 

CALL RANDU ( IX, IX , YFL) 

CGST=YF L"50. 

FC=(XI"CGST) 

2 S l M=0 . 

SSQ=0. 

DC 3 1=1, N 
SLM=SUM+DEM( I ) 

2 SSC = SSQ+DEM( I 2 
BAR=SUM/FLGAT (N) 

IF(BAR. LE.C.OJGO TO 16 
V AR= ( SSQ-6AR* SUM ) /FLOAT IN- l ) 

A = E AR- e * 2 / V A R 
ALPHA=B AR/ VAR 
DC 4 1=2,300 
Z( I J=2( 1-1 ) + l . 

XX=Z ( I - 1 ) -•‘ALPHA 

CALL G A M A ( A , X X , G A M , 6 , c R ) 

Y=1 .-GAM 
I F ( Y-PT ) 4 , 5, 5 

4 CONTINUE 

5 RN=Z( I- L ) 

IF(RN«LT.1.0) RN= 1 . C 
ALAMD=E>AR*12 . 

XNLM=2" ( AL AMC“ AX ) 

DNGM= X I *C 0 ST 

CVy = SQRT (XNUM/DNOM ) 

XLAM=ALAM0=5 . 

I F ( Q W • G T .X LAM ) C W = X L A M 
D EAR=AL AMC/4 . 

G=AMAX1( 1. C, ON, DEAR) 

• IF (K .LT .N) GO TO 7 
I AN=RN+YFL" ( C + L . ) 

7 M=N+ 1 

I F ( IAN-CEM CM > ) 8 , 9, 5 

8 SHRT= ( D EM ( M ) - 1 AN ) 

TSHT=TSFT+SHRT 
ACF= FLO AT ( IAN )/2. 

HCT=HCT+ACH=FC 

I ELY=RN/Q 

EUY = FLOAT ( IBUY + 1 ) *Q 
I AN=BU> 

GG TO 16 

9 CF= IAN-CEM ( M ) 

A C H= ( I AN + OH) /?. . 

FCT=HCT+AOH-HC 
IFCCH.GE.RNJGCj to 15 
IBLY=( RN-OHJ /Q 
BUY = FLOAT ( IBUY + 1 )*Q 
I AN--&UY +Ch 
GO TO 16 
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15 I £ N =0H 

16 N = N+ 1 

TSUM=TSUM+DEM(M) 

IF(N-57)2,17, 17 

17 TFCT=T hCT+HCT 
XSUN‘=XSUM + TSUM 
XSFT=X$HT+TSHT 
6G TO 2 C 

5 C CUE=1 .-XSHT/XSUM 

WRITE! 6,200 )XSHT,THCT, CUE 

2 0 C FORMAT ( 10X , « X SHT = ' , F 9 . 2/ 10X , ' THCT = • , F 10 . 2 / 1 CX , 
*10X,» GU E= ' ,F6.4) 

STOP 

END 



l 
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FORTRAN PROGRAM FOR NORMAL (SMO, SD) 

DIMENSION IDEM! 57) ,DEM( 57) ,Z(30C) 

AX=25 . 

X S UM= 0 . 

X Sl-T=0 . 

THCT=0 . 

XI=.01 

2C READI5 , 10 , END=50 ) ( IDEM C I ) , 1 = 1 » 57 ) 

1C FORMAT (1215) 

D'C 500 1=1,57 
5C0 CEM( I ) = I D E M ( I) 

P T= . 5 
Z( 1 ) = 0. 

N =2 1 
M = N 
K = 2 1 
F CT =0 . 

T SHT = 0 . 

T S UM=0 . 

DEV=0 . 

I X= 1237 56 1 

CALL RANDU ( IX, IX, YFL ) 

CCST=YFL*5C. 

HC= ( XI ~COST ) 

2 S U M = 0 . 

S SQ = 0 . 

DC 3 1 = 1, N 
SLM=SUM+DEM( I ) 

2 S S C = S S Q +0 E M ( I ) = * 2 
6 AR=SUM/FL CAT { N ) 

I F ( BAR . LE.O.OJGQ TC 16 
VAR=( S S C-6 AR :i SUM ) / FLOAT ( N-l ) 

S=SGRT ( VAR ) 

R K = to ;x S + 8 A R 
IFIRN.LT. L.C)RN=1.0 
AlAMD=i}AR* 12 . 

XMM=2- (ALAMD-AX ) 

C N C M = X I ^ C 0 S T 
Cto=SQRT (XNUM/DNOM ) 

X LAM=AL AMD~’ ; 3 . 

I F(QW.GT.XLAM) QW = XLAM 
D EAR= ALAMD/4 . 

Q = A M A X 1 ( 1 . 0 , Q W » 0 6 A R ) 

IFIK.LT .N ) GO TO 7 
I AN = RN+YFL^(i+l . ) 

7 M=N+ 1 

I F ( IAN-CEM (M) ) 3, 9, 9 

8 SHRT=(DEM(M)-IAN) 

TSHT = T ShT+SHRT 

A CF= FLO AT ( IAN ) / 2 . 

HCT=HCT + AOFI ; ‘HC 
I fit Y=RN /Q 

' BUY=FLO AT ( IBUY+1 )*Q 
I A (\ = B U Y 
GC TO 16 

9 CF = I AN-DEiM ( M ) 

A C F = ( lAN + OFi) /2. 

FCT=HCT+AOH-HC 
IF(CH .GE.RN) GO TO 15 
I BLY= ( RN-GH ) /Q 
6LY = FL0 AT ( IdUY + l ) ’"-Q 
1 AN=3UY + CH 
GC TO 16 

15 IAN=OH 

16 N = N + 1 

T SLM=T SLM + DEM ( M ) 

I F C N— 57 ) 2 , 17, 17 

17 TF.CT = T|-CT+FCT 

X SLM= X SL3 + T SUM 
XSFT=XShT+TSHT 
GC TO 2 C 

5 C OLE= l .-XSHT/XSUM 
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NR 1TE( 6,200 )XSHT,THCT, CUE 

200 FORMAT ( 10X,'XSHT=',F9.2/10X,'THCT=',F10.2/1CX, 
*1CX, • OUE= • ,F6 .4 ) 

STOP 

END 
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FORTRAN PROGRAM FOR NORMAL (x,S) 

DIMENSION I DEM (57 ) , OEM (57) ,Z(300) 

A — C *3 23 
B=C.717 
AX = 25 . 

XSLM=0 . 

XSFT=0. 

T F CT=0 . 

XI=.01 

2 C READ(5, 1Q,END=50) ( I DE M ( I ) , I = 1 , 5 7 ) 

1C FORMAT (1215) 

DC 500 1=1,57 
5 OC CEM ( I }= IDEM C I) 

PT = . 5 
Z( 1 ) =0 . 

N = 2 1 
M = N 
K = 2 1 
HCT=0. 

TS FT=0 . 

TSLM=0. 

C EV=0 . 

IX=1237561 

CALL RANDU( IX, IXtYFL) 

CCST = YFL- t 50. 

HC=(XI' ; CCST) 

2 SLM=0 . 

SSC=0. 

DC 3 1 = 1, N 
SLM=SUM+DEM ( I ) 

3 SSC= SSO+DEM ( I )'^2 
B A R = S U M / F L 0 A T ( N ) 

IF(8AR»LE.O.O)GO TO 16 

V AR=( SS C-B AP*SUM) / FLOAT (N-l ) 
XY=A+B*(ALOG( BAR ) ) 

AM AD=EXP (XY J 
S = 1 .2 5-'* AMAD 
SMC=3AR 
60 TO 13 

12 CEV = A3S (DEM(M)-SMC ) 

AMAD=G . 2 : '-DEV + C.8*AMAD 
S=1 ,25«AMAC 

S M C = 0 . 2 ,; D E M ( IT ) + 0 . 8 ~ S M 0 

13 RN = N"S + SMO 
IFtRN.LT.l .0 ) RN= 1.0 
ALAMO=SMO*12. 

X M U M= 2 * ( A L A iT 0" A X ) 

DNCM=X I *0 0 ST 
QW=SQRT ( XNLM/DNOM ) 

XLAM=ALAMC : -3 . 

I F (CW.GT.XLAM) QW = XLAM 
DEAR=AL AMD/A . 

C = AM AX 1 ( 1 .0, GW, DbAR ) 

IF(K.LT . N ) GO TO 7 
I AN=RN + YFL“(Q+l. ) 

7 M = N + 1 

1 F( IAN-DEM (Ml )8 ,9 ,9 

8 SFRT= ( 0 EM ( M ) - I AN ) 

TShT=TSFT+SHRT 
ACH=FLCAT ( IANI/2. 
hCT=HCT+AOH«HC 

I EUY = RN/G 

BLY = FLCAT ( I8UY+1)*G 
I AN = BUY 
GO TO 16 
5 Qh=I AN-UEM ( M ) 

A C h= ( IAMCH/2. 

HC7=HCT +ACh”FC 
IF(OH.GE.RN)GO TO 15 
I EOY= ( PN-OH ) /Q 
BIY=FLQAT( IBUY+1 )*G 
I AN = BUY +0H 
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GC TO 16 

15 I AN=0 h 

16 K=N+1 

T$LM=TSGK+DEM( M ) 

I F ( N-57 ) 12, 17, L7 

17 TPCT=T ECT + HCT 

X S LM=X SLM+ T SUM 
XSHT=XSHT+TSHT 
G C TO 2 0 

5C OLE=L.-XSHT/XSUM 

WRITE{6,20C)XSHT, TFiCTfOUE 

2 00 FORMAT ( 1CX, * XSHT = ' , F9 . 2/ 10X , ' THCT = ' , F 10 . 2/ 1 CX , 
* 1CX, ' OLE= ' » F 6 . 4 ) 

STOP 

EKC 
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APPENDIX D 
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Item # 6 (AF #174) 
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EMPIRICAL CDF FOR ITEM #@ (AF #43) 
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Demand 



APPENDIX E 



FORTRAN PROGRAM AND RESULTS OF SIMULATION MODEL 



DIMENSION ISHIP(7),A17),SEA(7),RSEA(7),NFCRT(7), 
1MCP(71,NSC(7),NCUP(7),INC(7),IRK(7) 

I C AY = 0 
I X = 7 73 
R C AY = 7 2 C 
R = 7 

N I CPM=0 
RSCP=0 
R X E AR I = C 
R X SGR I = C 
RXEARN=C 
RXSGRN=C 
RF=RDAY/3C 
CC 10 1=1,7 
ISHP( I 1 = 0 
R S E A ( I ) =1 
R FCRT (I ) = 0 
IRC ( I ) = C 
I R K ( I ) = C 
R ICP ( I ) =C 
RSC(I)=C 
RCLF ( I ) =0 
1C CCRTINUE 

PEAC(5,lGCC)A(L),A(2),A(3),A(4),A(5),A(t>),A(7) 

10CC FCRRATt 7( F3.C ) ) 

R E AC ( 5 » 1500 ) SEA ( 1) ? S E A ( 2 ) » S E A ( 3 } » S EA ( 4 } , S E A ( 5 ) » S E A ( 6 } » 
1SEA (7) 

15CC FCFRAT ( 7 (F5 .0 ) ) 

2 5 C CCRTINUE 

CC 15 J = l » 30 
1CAY=ICAY+1 

IF ( IDAY .GT .RCAY) GC TO 500 
CC 20 I = 1 » R 

IFtICAY.LT .NSEA ( I ) ) G 0 TC 18 
CALL RARDU(IX,IY,YFL) 

I X = I Y 

IRC ( I ) = - S E A ( I I-ALCG ( YFL ) 

R FCRT (I ) = ICAY + IR!C(I) 

R S E A ( I ) =NFCRT ( I ) + 1 C 

16 I F ( ISFIPt I I.GT.IDAY.CR.IShlP(I) .LE.OIGO TC ICC 
R ICP ( I ) = N I C P ( I ) + 1 
R CLP ( I)=NCLMt I) + l 

IF(ISHIF(I).GE.NPCPT(I).AND.ISHIF(I).LE.RSEA(I)), 

*GC TO 5 C 

CALL RAROUt IX, I Y, YFL) 

I X = I Y 

IRK ( I )=-A (I )«ALOG ( YFL ) 

ISFIPt I ) = I SI- IP ( I ) + IRK (I) 

GC TO 2C 

5 C RSC ( I ) = RSC ( I ) +NCUR ( I ) 

RCLP ( I ) =0 

CALL RARCL(IX,IY»YFL) 

I X = I Y 

IRK (I ) = -A ( I ) -ALCG ( YFL ) 

ISFIPt I ) = I S h I P ( I ) + I RK t I ) 

GC TC 2 C 

ICC IF ( ISHIFt I ) .GT. ICAY)GC TO 20 
CALL R A RCU ( I X f IY,YFL) 

I X = I Y 

IRK t I ) =- A t I )*ALOG( YFL ) 

ISFIPtI ) = ISt-IP( IMIRKtl) 

2 C CCRTINUE 
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15 CCMINUE 

RICP^=RICP< i)+RIICP(2)+MCP(3)+RI ICF(4)+N I CP <5HNICP(6), 
*4MCP(7 ) 

RSCN=RISC<1)+RISC(2)+RSC(3)+RISC<4MRSC(5)+RSC 16)+R!SC(7) 
kFITE(6 ,21CC 1MCPR ,RSCP 

2 1 C C FCPMAT< ‘0* , • ICP D E N AN D • , I 5 , 1CX , ' R S C DEMARC',15) 
RXEARI=RX6ARI+NICPP 
R > SCR I = RX SCR I +N I CPR*N I CPN 
NXEARN = RXBARN+NSClM 
R>SCRN = RXSCRN + RISCV*RSCf' 

CC 20 1=1,7 
R 1CP( I )=C 
RSC ( I ) =C 
2 C CCMINUE 
R ICFI* = C 
RSCF=0 
G C TO 2 c C 

5CC RXEARI=RXEARI/NP 
RXEARN=RXEARR/NP 
RXSCRI=NXSCRI/NP 
R>SGRN=RX$£RR/NP 
RVARI =NX S CRl-NXbAR I*NX EAR I 
RVARN=NXS£RN-NX6ARR»NXEARN 
PFITE(6 ,22 CO ) RXBAR I ,NVAPI,NXBARR ,RVARN 
2 2 C C FCPMAT ( '0 • , • ICP R E AN • , I 5 , 5X , 1 1 C P V AR I ANC E 1 , 15 , 

1 5 X , * NS C NEARI • ,I5,5X, 'NSC VAR I ANC E ' , 1 5 ) 

STCP 

ERC 



Variable names used in the program are defined as follows: 

ISHIP(I) Day of demand for Ship I 

NPORT(I) Day of arrival in port for Ship I 

NSEA(I) Day of departure from port for Ship I 

NICPM Total monthly demand recorded at Inventory 

Control Point (ICP) 

NSCM Total monthly demand recorded at Naval Supply 

Center (NSC) 

NXBARI Mean demand recorded at ICP 

NVARI ■ Variance of demand recorded at ICP 

NXBARN Mean demand recorded at NSC 

NVAFtN Variance of demand recorded at NSC 
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